function val = parzen(ds,Y,h)
     [m,n] = size(ds);
    val = [];
    for i=1:size(Y,1)
        p = Y(i,:);
        d = ds - repmat(p,[m,1]);
        d2 = sum(d.^2,2);
        val = [val; 1/size(ds,1)*sum(((h^2)*2*pi)^(-n/2)*exp(-((d2))/2/h))];
    end
end